
I Introduction 


The past decade has been marked by an increasing interest in performing 
active experiments in space* These experiments involve the artificial injec- 
ions of beams, chemicals, or waves into the space environment. Properly 
diagnosed, these experiments can be used to validate our understanding of 
plasma processes, in the absence 'of wall effects. Sometimes they even 
lead to practical results. For example, the plasma-beam device on SCATtlA became 
the prototype of an automatic device now available for controlling spacecraft 
charging at geostationary orbit. 

In this paper we discuss the future possibility of actively testing our 
current understanding of how energetic particles may be accelerated in space 
or dumped from the radiation belts using intense electromagnetic energy from 
ground based antennas. The ground source of radiation is merely a convenience. 
A space station source for radiation that does not have to pass through the 
atmosphere and lower ionosphere , is an attractive alternative. The text is 
divided into two main sections addressing the possibilities of (i) accelerat- 
ing electrons to fill selected flux tubes above the Kennel-Pet scheck limit 
for stably trapped fluxes and (2) using an Alfven maser to cause rapid deple- 
tion of energetic protons or electrons from the radiation belts. Particle 
acceleration by electrostatic waves have received a great deal of 
attention over the last few years ( Wong et al . , 1981; Katsouleas and 
Dawson , 1983 ). However, much less is known about acceleration using 
electromagnetic waves. The work described herein is still in evolution. 

We only justify its presentation at this symposium based on the novelty of 
the ideas in the context of space plasma physics and the excitement they have 
generated among several groups as major new directions for research in 
the remaining years of this century. 
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II Electron Acceleration by Electromagnetic Waves 

One of the first things we were mistaught in undergraduate physics is that 

electromagnetic (em) waves can't accelerate charged particles. If the particle 
gains energy in the first half cycle, it loses it in the second half. Teachers 
are, of course, clever people who want graduate students. So they hold off discus- 
sing gyroresonance, in which case, all bets are off. The resonance condition is: 


(1) w - k z v z - nfi 0 /y =0 


Here w is the frequency of the driving wave, k z the component of the wave vector 
along the zero order magnetic field _B 0 = _S 0 “z , v z the particle's component of 
velocity along Bo and n is an integer representing an harmonic of the gyro- 
frequency 2 o = q B Q /m, Y is the relativistic correction (1 - v^/c^) - !/^, 
q is the charge, and m the rest mass of the electron. 

Before going into a detailed mathematical analysis it is obvious that there 
are going to be problems accelerating cold ionospheric electrons to high energies. 
Higher than first gyroharmonics will have Bessel function multipliers where the 
argument of the Bessel function is the perpendicular covaponent of the wave 
vector and the gyroradius. For cold electrons with small gyroradii, all but the 
zero index Bessel function terms will be small. The second concern can be 
understood by considering the motion of a charged particle in a circularly 
polarized wave. Roberts and Buchsbaum (1964) have shown that with an electron 
in gyroresonance according to eq.(l) and _v j_ initially antiparallel to the wave 
electric field E and perpendicular to the wave magnetic field B, two effects 
combine to drive it away from resonance. As the electric field accelerates the 
electron, y increases, changing the gyrof requency . The magnetic component of 
the wave changes v z and thus, the Doppler shift term. It is only in the case of 
the index of refraction n = ck/ oj = 1 that unrestricted acceleration occurs. 

In all other cases the electron goes through cycles gaining and losing kinetic 
energy. 

Recently, the SAIC group ( Menyuk et al . 1986) has devised # a conceptually 
simple way to understand acceleration by era waves as a stochastic process. 

In terms of the relativistic moments p z and p j_ , eq.(l) can be rewritten as 

9 9 2 

P -(h - 1 ) p z + 2 r] z P z mc (n 2 J a) ) + ( (n Q J a> ) 2 - 1 1 me 2 

Depending on the phase velocity of the waves, equation (2) represents a family 
of ellipses ( n z = ck 2 / a) <1), hyperbolae ( n z >1 ) and parabolae 
( z = 1) in a p i , p z phase space. The zero order Hamiltonian can also be 
written in the form 

(2) H 0 /mc 2 = [ 1 + (p z /mc) 2 + (p j_ /me) 2 ] lj/2 - (p z /mc) ( to /ck z ) 
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Thus, in p j_ , p z space constant Hamiltonian surfaces represent families of 
hyperbolae ( n z < 1 ) ellipses ( n z > 1 ) and parabolae ( n z = 1 ). 
Hamiltonian surfaces have open topologies for indices of refraction n z < 1. 
The case n z = 1 in which resonance and Hamiltonian surfaces are overlying 
parabolae is that of unlimited acceleration studied by Roberts and Buschbaum 
(1964). 


In the case of small amplitude waves the intersections of resonance and 
Hamiltonian surfaces in p j_ , p z space are very sharp. As the amplitudes of 
ttie waves grow so too do the widths of resonance. For sufficiently large 
amplitudes, resonance widths may extend down to low kinetic energies allowing 
cold electrons to be stochastically accelerated to relativistic energies. 

It should be pointed out that although this model heuri sti cally explains 
the main conceptual reasons for stochastic acceleration to occur, its validity 
extends only to small angles 0 between k and _B Q . At large angles, it is 
not clear that the zero-order Hamiltonian topologies described above will still 
hold. 


Over the past several months we have developed a rigorous extension 
of the analytical model of Roberts and Buchsbaum by letting k. = k x ^ + k z £ 
assume an arbitrary angle to B > 0 . We begin with the Lorentz "equation. 

(3) d£ = q [ jE + _v x (_B 0 + _B) ] 
dt 

Tiie relativistic momentum and Hamiltonian are given by £ = n y v and 
H = mc^ y, respectively. The magnetic field of the wave J3 i s related to the 
electric E through Maxwell's equation B_ = (c/w )k. x _E. The time rate of change 
of the Hamiltonian is 

• 

(4) H = q EM v_ = qc^ E. p/H 

If we define E x = E^ cos <f> , E v = - E 2 sin b and E z = - E 3 cos <J> , where 
b - k x x + k 2 z - m t then equation (4) may be rewritten in the form 

. 

(5) H H = ^Ej_ p x cos ^ Py siri <> - p z cos ^ 

c z co oj a to 


The Lorentz force equation can also be rewritten as 


(6) Px +Py [ ^ + qE? sin (J) J = qE] ( to - K z z) cos $ 

my w oj 


(7) Py - Px [ ^ + qE? ^y sin 4> ] = _ qE? (u) _ k z z) sin $ 

my a) a) 

• • • 

(8) P z - K^_ H + E3_ (p x +n p ) = 0 

u E : 


where K z - k z (1 + E 3 k x /E 1 k z ). Equations (5-8) are exact. Our first simplifi- 
cation is to assume E 2 k x / u> = B z << B Q , then eqs. (6-8) may be combined 
to give 

* t 

(9) 4HH = q ( E^ + E 2 ) [ / Q T cos( a + $ - a T + $ f ) dt r + 

0) 03 o 

r t 

+ j R 1 cos ( a + $ - a f - $ f )dt f - 2p j_ sin ( a + j> + a ) ] 

0 

r f t 

+ q ( E]_ — E2 ) [ J Q 1 cos ( <p - o + o r - <p ') dt T 
0) o 

,t 

+ J R' co s( 4> - a + 6 '+ a ') dt’+ 2p j_ sin ( 6 - a - a ) J 
0 

- E 3 { 4 ( p zo + K 7 (H-II 0 ) ) COS 4) 

0) 03 

li / (Q f + R ! ) [ cos( $ + <{> ') + cos ( <t> ~ <t> f ) ] dt T 1 

Ei o 

where a ( t) = /* fl ( f) df, tan a = - ( p xo /p yo ) , 

(the subscript o refers to the initial conditions at t = 0), and 
Q = q E 1 ( m - K z z) - qE? ( w - k z z) 

0) 03 


• • 

R = qE] ( m - K z z) + qE? ( u - k z z) 

W tj) 

Primed and unprimed quantities are evaluated at times t' and t, respectively. 

We note that accelarations represented in Eq. (9) are related to terms multiplying 
electric fields in right-hand (E t + E 2 ), left-hand (E], - E ? ) and parallel E 3 model. 
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Our next simplification is to substitute for x and z in eq.(9) the zero 
order solutions (in the electric field amplitude) of eqs. (6-8), That is, we 
take x = p cos( o + a ) where p = v j_ / Q is the electron gyroradius and 


Ou) Pz = [ ,p zo + k_2_ (h - u 0 ) j . 

00 

We note that eq.(10) reduces to eq.(2) by taking K z = k z , which is only valid 
for small angles between _k and _B 0 . In fact, Figure 1 shows that Hamiltonians 
with open (hyperbolic or parabolic) topologies in p z , p j_ space at small 
angles between k. and become closed (elliptical) as the angle increases* 

The practical implication is that cases of potentially infinite acceleration 
with k - k 2 become restricted to finite values at other directions of wave 
propagation. 

By taking x = p cos( O + a ) and expanding terms with sin k x x and 
cos k x x in series of Bessel functions, eq . (9) becomes 


(11) 4HH = l T n 

c~ oo n 

T n = ^ E 1 + e 2) J n-1 ( k x P ) { 1 ! [ Q f J f cos(n B+mG'+^ + iJ;') 

oo mo m+1 

+ cos(n 0 - m 0 ’ + ip — if; ') ] dt' + 2p j_ cos (n 0 + ™ ) } 

+ jL_ C 1^1 — E 2 ) Jn+1 (k x P ) { \ j [ Q'J' cos(n 0 - m 0 '+ - i(; ' ) 

oo m o nrt 1 

+ k! J f ni-l cos ( n $ + m 6 f + ^+^ f )j dt f + 2p j_ cos (n 9 + ) J 

“ p ) { ^ ( P Z o + ) cos (n 0 + ) 

w CJ 

- E, l j 1 (Q' + R ' )J ' [ cos( n 6 + m 0 '+ i|» + i|/ * ) 

p m o m 


+ cos(n 0 -mB'+^-iJ; 1 )] dt f } 


where 0 = / t fi(c')dt , +a+Tr/2, J' = J (k p '), ( v = m, ra + 1) 

o v v x - 

and = k z z - oo t. 

After averaging over the fast (gyroperiod) time dependencies and a good deal 
of tedious algebra, we obtain that, for each n, the particle energy obeys the 
following differential equation: 

(12) (U + l) 2 ( 1_ dU ) 2 + V n (U) = 0 

w dt 

where U = (H-H 0 )/H 0 and 
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2 


2 

v n ( ,J ) = U 2 ( U +2 r n /d 1 
4 


^ (o) sin 4> n d 1 U ( U + 2r n /d l ) 


+ 2 1 ~ l 2 { ( 1 2 d l“ z l h l) ( G n+1 ( u ) + F n+l(U) ) 
2 

+( E 2 d 2 - E x h 2 ) F n+1 (U) } 


- E i + E £ { ( E l d i + ^ 2 d l^ ( G n-l( d ) + Fn-i(U) ) 

2 

+ ( E 1 h 2 + E 2 d 2 ) f n _i(U) } 

2 2 
- E 3 { h L ( G n (U) + F n (U) ) + h 2 F n (U) } - ( ♦ (O)cos $ n ) 

where E i = - (q E^/ u ) c/H Q (i=l ,2,3) , d^ = 1 - K z k z c^/ to ^ 

d 2 = K z k z c2/ w 2 " k z z o/ w » h l " 1 + K z /k z (d l " O 

r n = 1 ” k z z o/ “ “ n « Q / to , h 2 = K z /k z d 2 

'f' (0) = v x q /2c [ ~( E 3 + E 2 ) J n _ L (k x P 0 ) + ( 2 2 ~ 2 i> 

d n+l (k x P o) ] + v z0 /c E 3 J n (k x p Q ) , 

4 > n = n ( a +_jr_) + k z z Q 
2 

,U 2 r i 

and G v (U) = J J v [ k x p (U 1 ) J U’ dU’ 

o 

,U 2 f , 

F v (U) = / J v [ k x p CJ’) j dU’ , ( v = n, n + 1). 
o 

Eq.(l'2) is in the form of the equations of a harmonic oscillator. Under the 
limit 0=0, Eq. (12) becoiaes the equation derived by Robert and Buchsbaum 
(1964). The limits of the particles excursion in energy for a given resonance 
n and electric field E can be found by setting the potentials V n (U) = 0. At 
wave amplitudes where the range of potentials for different harmonics overlap, 
we have the onset of stochast icity . 

At the present time we have just begun to explore the numerical solutions of 
equation (12). In Figure 2, we show some of our preliminary results. We assume 
that w p e / ft 0 = 0.3, the electric field amplitude is such that Si = 0.1, and the 
wave frequency is a) = 1.8 ft 0 . We consider only the second cyclotron harmonic 
since this is the closest to satisfying the resonance condition, eq.(l), for 
initially cold electrons. The components of the wave electric field and the 
refractive index ri are calculated from the cold plasma dispersion relation for 
electromagnetic waves at any arbitrary angle 0 to B Q . It turns out that rt is 
always smaller than, but very close to 1 ( n = 0.97). The maximum allowed 
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energy gain, as given by the zero order Hamiltonian topologies, is represented 
by the solid lines • The shaded region represents the actual energy gain as 
obtained by requiring (U) < 0. We see that for 0 - 35°, initially cold 
electrons can be accelerated to very high energies. In fact, for cold 
electrons we find that U = Y - 1 and that the particle can gain as much as 2.3 
Mev. As 0 decreases more initial kinetic energy is required for any 
acceleration to take place. For large 0 , the elliptical haniltonian 
topologies severely restrict the energy gain. 


Ill The Alfven Maser 


Active control of energetic particle fluxes in the radiation belts has 
maintained a continuing interest in both the United States and the Soviet Union. 
Electron dumping experiments concluded by the Stanford University and Lockheed 
groups using VLF transmissions are well known ( Inan et al. 1982, Imhof et al . 
1983). Perhaps less known is a theoretical paper by Tr akthengert s (1983) 
entitled "Alfven Masers" in which he proposes a theoretical scheme for dumping 
both electrons and protons from the belts. The basic idea is to use RF energy 
to heat the ionosphere at the foot of a flux tube to raise the height integrated 
conductivity. The conductivity is then modulated at VLF or ELF frequencies 
which modulates the reflection of waves that cause pitch angle diffusion in the 
equatorial plane. The artifically enhanced conductivity of the ionosphere thus 
maintains high wave energy densities in the associated flux tube, thereby 
producing a masing effect. 

In addition to external ionospheric perturbations particle precipitation 
also raises ionospheric conductivity. The masing of the VLF waves causes 
further precipitation which, in principle, results in an explosive instability. 
The purpose of this section is to establish the basic equations and to present 
the results of a preliminary computer simulation. 

The fundamental equations derived by Tr akhtengert s (1983) are based on 
quasilinear theory and relate only to the weak diffusion regime. It is useful 
to use a similar set of equations derived by Schulz (1974) based on phenomeno- 
logical arguments that includes strong pitch angle diffusion. The key variables 
are N, the number of trapped particles per unit area on a flux tube and £ the 
wave intensity averaged over the flux tube. In this we assume that e is 
directly proportioned to the pitch angle diffusion coefficient. The time rate 
of change for N is 

(13) dN —A £ N + g . 

dt 1 + e t 

where the first term represents losses due to pitch angle scattering with A a 
constant and S represents particle source terms in the magnetospheric 
equatorial plane, t is a parameter that characterizes lifetimes against 
strong pitch angle diffusion. The time rate of change of e is given by 

(14) d £ = [ 2 y *N/N* 1 £ + Vg £ In R + W 

dt 1 + ex LRg 


220 



ORIGINAL PAGE IS 
OF POOR QUALITY 


The first term represents wave growth near the equatorial plane, the second term 
gives the wave losses in and through the ionosphere and the third accounts for any 
wave energy sources. The terms Y * and N* are used to denote the weak diffusion 
growth rate and column density of a flux tube at the Kennel and Petschek (1966) 
limit for stably trapped particles. In the second term, v g /LR e approximates 
bounce frequency of waves where v g is the group velocity of the wave LR the 
approximate length of a flux tube; R is the reflection coefficient of the ionos- 
phere. Since R < 1 the second term is always negative. The (1 + e x ) term 
empirically lowers growth rate due to the pitch angle distribution becoming more 
isotropic under strong diffusion conditions. 

In our present study we have examined numerical solutions of equations 
(13) and (14) using non-equilibrium initial conditions. The first case is 
represented by Figure 3 in which we started initial wave energy densities 
wnich are a factor of 3 (top panel) and 0.1 (bottom panel) above the Kennel— 
Petschek limit. In both cases we ignored associated enhancements in ionos- 
pheric coupling that lead to increased reflectivity. We see that the wave 
density quickly damps to the Kennel— Petschek equilibrium represented 
by the solid line. 


In the second level of simulation the wave energy density is initially set 
at a factor of three above the Kennel— Petschek equilibrium value but includes 
a coupling factor to the ionosphere 4 . We find that for values of 4 > 10% 
the oscillations become spike— like. The top panel of Figure 4 represents the 
normalized wave energy density for 4 = 10% after the waves have evolved into 
periodic spitces. ihe middle and bottom panels of Figure 4 represent the nor- 
malized energetic particle density (cm"2) contained on a flux tube and the 
normalized height integrated density of the ionosphere. Attention is directed 
to the phase relationship between the maxima of the three curves. The maxi- 

mum. pnprcrpfi'n nprfihlo flnv 1 a nA c tbo ~ , — — a — *-i i, ^ i 

* o r — ^ ^ ^ ^ x v- < j.v* o cut. wave tcLm ai iu gucb LiiLuugii l lie tvenne-L~ 

Petschek value as the wave growth changes from positive to negative. 



Fig. 3. Example of wave energy densities initially set at factors 
of 3.0 and 0.1 above Kennel Petschek equilibrium value. 
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Fig. 4. Example of spike-like wave structures as well as energetic 
particle losses and ionospheric density changes with 

magnetosphere-ionosphere coupling . 
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Fig. 5. Simulated, raormalized wave energy density with magnetosphere- 
ionosphere coupling. A VLF source is turned on at t = 650s. 
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The maximum ionospheric effect occurs after the wave spike maximum. Our phys- 
ical interpretation of Figure 4 is as follows. A spike in the wave energy 
density causes a depletion of electrons trapped in the belts to levels well 
below the Kennel-Petschek limit. The subsequent drop of precipitating electron 
flux allows the ionospheric conductivity to decrease. Thus, VLF waves are less 
strongly reflected back into the magnetosphere. This effectively raises the 
Kennel-Petschek limit as higher particle fluxes are necessary to offset in- 
creased ionospheric VLF absorbtion. In the presence of equatorial sources of 
particles, the similations show flux levels building to 1.15 times the Kennel- 
Petschek limit. The enhanced fluxes in the magnetosphere, even with weak pitch 
angle diffusion, allows the ionospheric conductivity to rise, eventually leading 
to another masing spike. 

Figure (5) shows the effect of an external VLF signal. The first few spikes 
result i- rom the masing effect of the ionosphere due to particle precipitation. 

At t = 650 seconds a VLF square wave source is turned on with a 50 second dura- 
tion. The spikes now are modulated at the driving frequency at a reduced ampli- 
tude. The amplitude is reduced since the fluxes are more frequently dumped 
with the VLF signal present than in its absence. 

Ivers en et al . (1984) using simultaneous ground and satellite measurements, 
have recently observed the modulation of precipitating electrons at pulsation 
frequencies. In terms of our simulations these would be close to the situation 
shown in Figure 4 in which natural masing occurs in a flux tube. The observed 
frequencies are consistent with those expected from the linear theory. Detailed 
comparison with experimental data necessitates knowing the efficiency with which 
VLF waves reach the ionosphere. 


IV Conclusion 


Although the work presented in this paper is still in a very preliminary 
stage of development it appears that significant space effects can be produced 
by the injection of intense electromagnetic waves into ionospheric plasmas. 

In the coming months we expect that as calculations mature we will grow in the 
ability to translate mathematical representation into physical understanding. 

If the results of our analyses live up to early promise then a series of ground- 
based wave emission experiments will be developed to measure injection effects 
in space. The upcoming ECHO— 7 experiment presents a well instrumented target of 
opportunity for electron acceleration experiments with the HIPAS system. After 
the launch of the CRRES satellite it will be possible to make simultaneous in 
measurements of wave and particle fluxes in artificially excited Alfve~ 
Masers. Looking forward to the 1990 f s it appears that the WISP experiment 
planned for the Space Station will make an ideal source for both electron 
acceleration and radiation belt depletion experiments. Recently a Soviet 
experiment measured electrons accelerated to kilovolt energies using a low 
power telemetry system ( Babaev et al . , 1983). Just image that what could be 
done with the specifically designed, high power WISP I 
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